PL-TR-94-2243 


REMOTE  SOUNDING  OF  CIRRUS  CLOUD  PARAMETERS 
DURING  DAYTIME 


S.C.  Ou 
K.N.  Liou 
N.X.  Rao 


Center  for  Atmospheric  and 
Remote  Sounding  Studies  (CARSS) 
Department  of  Meteorology 
University  of  Utah 
Salt  Lake  City,  Utah  84112 


1  September  1994 


Final  Report 
1  May  1992-1  June  1994 


LO 


Approved  for  public  release;  distribution  unlimited 


LO 


PHILLIPS  LABORATORY 
DIRECTORATE  OF  GEOPHYSICS 
AIR  FORCE  MATERIEL  COMMAND 
HANSCOM  AFB,  MA  01731-3010 


DTI*  QUAUTT  IKSParaSD  B 


This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


VINCENT/  FALCONE 
Contract  Manager 


DONALD  A.  CHISHOLM,  Chief 
Satellite  Analysis  and  Weather  Prediction  Branch 
Atmospheric  Sciences  Division 


Director,  Atmospheric  Sciences  Division 


This  document  has  been  reviewed  by  the  ESC  Public  Affairs  Office  (PA)  and  is  releasable  to  the 
National  Technical  Information  Service  (NTIS). 


Qualified  requestors  may  obtain  additional  copies  from  the  Defense  Technical  Information  Center 
(DTIC).  All  others  should  apply  to  the  NTIS. 


If  your  address  has  changed,  if  you  wish  to  be  removed  from  the  mailing  list,  or  if  the  addressee 
is  no  longer  employed  by  your  organization,  please  notify  PL/m  ,  29  Randolph  Road,  Hanscom 
AFB,  MA  01731-3010.  This  will  assist  us  in  maintaining  a  current  mailing  list. 


Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notices  on  a  specific 
document  require  that  it  be  returned. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


Public  reporting  burden  for  this  collealon  of  information  is  estimated  to  average  i  hour  per  response,  including  the  lime  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  colleaion  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspea  of  this 
collection  of  information,  including  suggestions  for  r^ucing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  information  Operations  and  Reports.  1215  Jefferson 
Oavis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  and  to  the  Office  of  Management  and  Budget.  Paperwork  Reduction  Project  (0704-0 188),  Washington,  DC  20503. 


1.  AGENCY  USE  ONLY  (Leave  blank) 


4.  TITLE  AND  SUBTITLE 


2.  REPORT  DATE 

1  September  1994 


3.  REPORT  TYPE  AND  OATES  COVERED 

Final  Report/  1  May  1992  -  1  June  1994 


5.  FUNDING  NUMBERS 


Remote  Sounding  of  Cirrus  Cloud  Parameters  during  Daytime 


6.  AUTHOR(S) 

S.C.  Ou/  K.N.  Liou/  and  N.X.  Rao 


F  19628-92-K-0019 
PE  61102F 

PR  2310  TR  CP  WU  KA 


7.  PERFORMING  ORGANIZATION  NAME(S}  AND  ADORESS(ES) 

Center  for  Atmospheric  and  Remote  Sounding  Studies  (CARSS) 
Department  of  Meteorology 

University  of  Utah/  Salt  Lake  City/  Utah  84112 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  ADORESS(ES) 

Phillips  Laboratory 

29  Randolph  Road 

Hanscom  AFB/  MA  01731-3010 

Contract  Mainager:  Vincent  Falcone/GPAB 


10.  SPONSORING /MONITORING 
AGENCY  REPORT  NUMBER 


PL-TR-94-2243 


12a.  DISTRIBUTION /AVAILABILITY  STATEMENT 


12b.  DISTRIBUTION  CODE 


Approved  for  public  release;  distribution  unlimited 


13.  ABSTRACT  (Maximum  200  words) 


A  numerical  scheme  has  been  developed  to  remove  the  solar  component  in 
the  AVHRR  3.7  /urn  channel  for  the  retrieval  of  cirrus  parameters  during 
daytime.  The  cloud  retrieval  scheme  described  in  Ou  et  al.  (1993)  has  been 
modified  in  connection  with  the  removal  program.  We  have  applied  the 
removal -retrieval  scheme  to  the  AVHRR  daytime  data,  collected  during  FIRE-I- 
IFO.  Distributions  of  the  retrieved  cloud  heights  and  optical  depths  are 
comparable  to  those  determined  from  lidar  measurements  and  from  GOES  visible 
and  IR  channels  data  reported  by  Minnis  et  al.  (1990a).  Moreover,  validation 
of  the  algorithm  has  been  performed  by  using  various  data  sets  that  were 
collected  during  the  FIRE-Il-IFO  (November  -December,  1991)  at  Coffeyville, 
Kansas.  We  have  focused  on  the  5  December  case.  We  demonstrate  that  the 
retrieved  cirrus  cloud  temperature,  mean  effective  ice  crystal  size,  And 
optical  depth  closely  match  the  values  derived  from  the  balloon-borne 
replicator  measurements . 


14.  SUBJECT  TERMS 


Cirrus  clouds/  remote  sensing/  AVHRR 


15.  NUMBER  OF  PAGES 
68 


16.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION  18.  SECURITY  CLASSIFICATION  19.  SECURITY  CLASSIFICATION  20.  LIMITATION  OF  ABSTRACT 
OF  REPORT  OF  THIS  PAGE  OF  ABSTRACT 

Unclassified  Unclassified  Unclassified  Unlimited 


NSN  7540-01-280-5500 


Standard  Form  298  (Rev.  2-89) 

Prescribed  by  ANSI  Std  Z39-18 
298-102 


TABLE  OF  CONTENTS 


Page 


1.  INTRODUCTION . 1 

2.  REMOVAL-RETRIEVAL  SCHEME  . 5 

2.1  Detection  of  Clear  and  Cloudy  Pixels  . 5 

2.2  Removal  of  3.7  /im  Solar  Component  . 7 

2.3  Retrieval  of  Cirrus  Cloud  Parameters  . 15 

2.4  Removal -Retrieval  Scheme  . 18 

3.  SENSITIVITY  STUDIES  AND  ERROR  ANALYSES  . 21 

4.  APPLICATION  TO  SATELLITE  DATA  ASSOCIATED  WITH  FIRE-I-IFO  . 32 

4.1  Retrieval  over  a  l'^  x  1°  Area  . 34 

4.2  Retrievals  over  a  5°  x  5°  Area  . 39 

4.3  Retrievals  over  a  20°  x  30°  Area  . 44 


5.  VERIFICATION  USING  FIRE-II-IFO  MEASUREMENTS  . 50 

5.1  Data  Sources  . 50 

5.2  Results  and  Discussions  . 51 


6.  CONCLUSIONS 


REFERENCES 


iii 


.60. 


.62 


o  ^ 


;  B.y' _ 


i  1. 1 OD 


_ _ Codes 


:Mat 


:\ 


SpGoiei 


. .  "  - 


1. 


INTRODUCTION 


Cirrus  clouds  are  global  in  nature  and  primarily  exist  in  the  upper 
troposphere  and  lower  stratosphere.  These  clouds  are  composed  almost  entirely 
of  ice  crystals.  Cirrus  clouds  play  a  key  role  in  the  global  energy  balance 
through  latent  heat  and  radiative  energy  exchanges.  Information  on  cirrus  cloud 
parameters  is  critically  important  to  the  development  of  cirrus  forecast  models, 
and  the  upgrading  of  real-time  global  cloud  analysis  (Liou  1986). 

Inference  of  the  microphysical  and  optical  properties  of  cirrus  clouds 
using  satellite  or  aircraft  data  has  been  reported  by  a  number  of  researchers. 
Wielicki  et  al.  (1990)  used  the  Landsat  0.83,  1.65,  and  2.21  ^m  channel 
reflectances  and  phase  functions  for  ice  crystals  to  infer  cirrus  cloud  particle 
sizes,  assuming  that  single  scattering  dominates  the  radiative  transfer  within 
cirrus  clouds.  Liou  et  al.  (1990)  developed  a  physical  retrieval  method  to  infer 
the  temperature  and  optical  depth  of  tropical  cirrus  anvils  using  the  data  of  the 
dual  channel  (6.5  and  10.5  /im)  ,  downward -viewing,  narrow-field-of-view 
radiometers  on  board  NASA  ER-2.  Minnis  et  al .  (1993a,  b)  developed  a  technique 
based  on  the  radiative  transfer  theory  and  parameterizations  to  retrieve  cirrus 
optical  depth  and  height  using  the  Geostationary  Operational  Environmental 
Satellite  (GOES)  0.65  and  11.5  fim  channel  radiances.  Ou  et  al.  (1993)  presented 
a  physical  retrieval  scheme  using  the  Advanced  Very  High-Resolution  Radiometer 
(AVHRR)  3.7  and  10.9  /urn  channels  for  the  inference  of  nighttime  cirrus  cloud 
parameters,  including  cloud  temperature,  optical  depth,  emissivity,  and  mean 
effective  ice  crystal  size  based  on  the  theory  of  radiative  transfer  and 
parameterizations.  This  retrieval  scheme  has  been  applied  to  the  nighttime  AVHRR 
satellite  data  collected  at  0930  UTC,  28  October,  1986,  over  the  region  of  the 
First  ISCCP  Regional  Experiment,  Phase  I,  Cirrus  Intensive,  Field  Observation 
(FIRE-I-IFO) .  The  retrieved  cirrus  height  and  mean  effective  ice  crystal  size 
compare  reasonably  well  with  ground-based  lidar  and  aircraft  measurements. 
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During  daytime,  the  3,7  fim  channel  radiance  contains  both  solar  reflection 
and  thermal  emission.  To  apply  the  retrieval  scheme  to  the  daytime  satellite 
data,  the  3.7  /im  channel  solar  component  must  be  quantified  and  removed.  The 
problem  of  removing  the  3.7  /im  solar  reflection  (referred  to  as  "sol")  and 
thermal  infrared  emission  (referred  to  as  "IR")  component  for  remote  sensing 
applications  has  been  investigated  by  several  researchers  in  the  past.  Arking 
and  Childs  (1985)  included  the  solar  component  of  the  3.7  /im  channel  in  the 
formulation  of  a  general  retrieval  program  for  cloud  parameters.  Kleespies 
(1993)  proposed  a  dual -observation  method  to  remove  the  thermal  infrared 
component  in  the  3.7  /im  channel  radiances  for  inferring  the  particle  sizes  of 
opaque  marine  stratiform  clouds.  In  this  report,  we  describe  a  method  to 
remove  the  solar  component  in  the  AVHRR  3.7  /im  channel  radiance  for  the  retrieval 
of  cirrus  cloud  parameters  during  daytime.  This  method  involves  various 
threshold  values  for  AVHRR  Chs.  1  (0,63  /im)  ,  2  (0.8  /im)  ,  3  (3.7  /im)  ,  4  (10.9  /im)  , 
and  5  (12  /im) ,  which  are  prescribed  to  separate  clear  and  cloudy  pixels.  A  look¬ 
up  table  relating  the  Ch.  1  reflectance  to  the  Ch.  3  solar  reflectance  is 
subsequently  constructed.  This  table  is  based  on  the  effective  surface  albedos, 
selected  ice  crystal  sizes,  and  a  set  of  sun-satellite  geometric  parameters.  A 
comprehensive  adding -doubling  radiative  transfer  program  is  used  to  generate 
numerical  values  in  the  look-up  table.  Removal  of  the  Ch.  3  solar  component  is 
accomplished  by  using  the  look-up  table  and  the  measured  Ch.  1  reflectance.  For 
retrieval  of  cloud  parameters  during  daytime,  the  scheme  described  in  Ou  et  al. 
(1993)  is  modified  in  connection  with  the  removal  program. 

Validation  of  the  satellite  remote  sensing  technique  for  the  retrieval  of 
cirrus  cloud  parameters  requires  ground  and  air- truths  that  are  available 
from  a  composite  field  experiment.  In  addition  to  our  efforts  cited  above, 
Minnis  et  al .  (1990b;  1993b)  have  used  ground-based  lidar  measurements  collected 
during  FIRE-I-IFO  to  verify  the  retrieved  cirrus  cloud  height.  Baum  et  al. 


2 


(1994)  have  also  presented  comparison  results  for  the  satellite  retrieved  and 
lidar,  radar,  and  aircraft-measured  cloud  heights,  using  data  collected  during 
FIRE-II-IFO.  However,  attempt  to  validate  the  retrieved  optical  depths  and  ice 
crystal  sizes  using  aircraft  data  has  not  been  made  at  this  point.  The  prime 
difficulties  for  this  are  due  to  the  horizontal  nature  of  aircraft  flight  path 
and  the  inhomogeneity  of  cirrus  cloud  properties.  To  evaluate  the  cirrus  cloud 
vertical  optical  depth,  continuous  microphysical  measurements  in  the  vertical 
must  be  available.  Moreover,  to  minimize  the  effects  of  horizontal 
inhomogeneity,  the  cloudy  region  sampled  by  the  in-situ  instruments  must 
correspond  to  that  viewed  by  the  satellite. 

We  have  carried  out  a  more  comprehensive  validation  for  the  retrieved 
cirrus  cloud  optical  depths  and  mean  effective  ice  crystal  sizes,  using  the  NOAA- 
11  AVHRR  data,  the  balloon-borne  replicator  data,  and  the  NCAR-CLASS  (Cross-chain 
LORAN  Atmospheric  Sounding  System)  sounding  measurements  that  were  obtained 
during  FIRE-II-IFO.  The  advantage  of  using  the  replicator  data  is  that  the 
replicator  sondes  will  provide  a  continuous  record  of  the  cloud  microphysical 
properties  in  the  vertical.  Thus,  it  is  possible  to  determine  the  optical  depth 
and  vertically  averaged  mean  effective  ice  crystal  size.  Furthermore,  we  have 
also  used  the  retrieved  cirrus  cloud  optical  depths  and  mean  effective  ice 
crystal  sizes  in  conjunction  with  a  radiative  transfer  parameterization  program 
to  compute  the  surface  radiative  fluxes  which  are  subsequently  compared  with  the 
ground-based  radiometer  measurements  that  were  obtained  during  FIRE-II-IFO. 
These  comparisons  will  thus  provide  a  consistent  check  on  the  retrieved  cloud 
parameters,  on  the  one  hand,  and  the  radiation  parameterization  program,  on  the 
other. 

This  report  is  organized  as  follows.  Section  2  describes  the  clear/cloud 
detection  scheme  and  the  removal -retrieval  procedures.  Section  3  presents  the 
results  of  sensitivity  studies  and  discusses  the  effects  of  various  possible 
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error  sources  on  the  accuracy  of  removal -retrieval  results.  In  Section  4,  we 
demonstrate  the  applicability  of  the  retrieval  algorithm  to  satellite  data 
associated  with  FIRE-I-IFO.  Section  5  then  presents  results  of  a  validation 
effort  using  data  obtained  from  two  cases  of  FIRE-II-IFO.  Finally,  conclusions 
are  given  in  Section  6. 
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2. 


REMOVAL-RETRIEVAL  SCHEME 


2.1  Detection  of  Clear  and  Cloudy  Pixels 

The  clear/cloud  detection  scheme  is  based  on  several  criteria  described  in 
Saunders  and  Kriebel  (1988).  Figure  1  presents  a  schematic  flow  chart  of  the 
detection  scheme  for  identifying  clear  pixels  during  daytime.  The  first 

criterion  involves  the  use  of  the  Ch.  4  brightness  temperature  (T^)  .  For  a  clear 

pixel,  T^  must  be  greater  than  a  defined  threshold  temperature  (T^c)  .  We  may  use 

the  map  of  over  the  clear  portion  of  the  data  domain  to  derive  a  mean  clear 

brightness  temperature  (T^)  .  The  T^c  is  empirically  set  as  T^  -  2  K.  Pixels  that 
satisfy  the  first  criterion  subsequently  undergo  the  second  test,  which  requires 
that  the  Ch.  1  reflectance  (r^)  be  smaller  than  a  threshold  value  (r^c)  ,  i.e.  r^ 
<  r^c*  The  threshold  r^c  is  determined  based  on  the  identification  of  peaks 
corresponding  to  cloudy  and  clear  conditions  in  the  histogram  of  r^. 

Clouds  generally  have  high  reflectances,  whereas  most  types  of  ground  and 
sea  surfaces  (except  snow  and  ice)  have  low  reflectances  at  the  visible 
wavelengths.  Using  this  property,  a  visible  threshold  may  be  set  up  to 

differentiate  between  cloudy  and  clear  reflectances.  Because  the  clear 

reflectance  varies  with  the  surface  type,  the  visible  threshold  is  a  function  of 
surface  characteristics.  The  third  criterion  uses  the  ratio  of  Ch.  2  (0.8  /xm) 
to  Ch.  1  (0.63  fim)  reflectances  (Q  =  r2/ri)  .  This  ratio  has  a  unique  feature 
with  respect  to  surface  types.  For  clouds,  Q  is  closer  to  1  because  the  cloud 
reflectances  for  these  two  wavelengths  are  about  the  same.  Over  water  surfaces, 
Q  <  1  because  the  reflectance  of  the  0.63  /xm  wavelength  is  more  pronounced.  Over 
land  with  vegetation,  Q  >  1  because  the  reflectance  at  the  0.8  /xm  wavelength  is 
much  larger  than  that  at  the  0.63  /xm  wavelength.  We  define  a  threshold  value 
such  that  a  pixel  is  clear  if  Q  >  Q^.  We  have  set  a  uniform  value  of  1.6 
(corresponding  to  land) .  This  threshold  value  is  applied  to  all  the  present 
retrieval  analyses,  because  most  of  the  area  within  the  retrieval  domains  are 
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Daytime  Clear/Cloud  Detection  Scheme 


Fig  1 ,  A  flow  chart  for  the  daytime  clear/cloud  detection  scheme 


over  land.  The  only  significant  water  surfaces  involved  are  over  the  Great  Lake 
region.  However,  the  display  of  Ch.  4  brightness  temperatures  shows  that  this 
area  is  mostly  cloudy.  Pixels  over  that  area  are  identified  as  cloudy  by  other 
criteria  (see  Section  4.c).  Finally,  the  fourth  criterion  uses  the  difference 
between  Ch.  4  and  Ch.  5  brightness  temperatures  which  must  be  less  than  a 
threshold  value,  6  (-2  K)  .  If  a  pixel  satisfies  all  four  criteria,  it  is 
identified  as  clear.  Otherwise,  it  is  classified  as  cloudy.  These  criteria  are 
complementary  to  each  other.  For  example,  low- cloud  pixels  are  identified  as 
cloudy  because  their  Ch.  1  reflectances  are  higher  than  the  threshold  (criterion 
2),  even  though  their  brightness  temperature  difference  between  Chs.  4  and  5  is 
smaller  than  2K  (criterion  4).  Also,  pixels  associated  with  thin  cirrus  are 
identified  as  cloudy  as  long  as  their  Ch.  4  brightness  temperatures  are  lower 
than  the  threshold  value  (criterion  1),  even  though  their  Ch.  1  reflectances  and 
brightness  temperature  differences  are  lower  than  threshold  values. 

2.2  Removal  of  3.7  /im  Solar  Component 

Let  the  angle  of  the  incident  solar  radiation  be  Oq,  while  the  satellite 
viewing  angle  be  $,  The  relative  azimuthal  angle  (A<^)  is  defined  as  the  angle 
between  the  horizontal  projection  planes  of  the  sun  and  the  satellite.  We  define 
fiQ  and  to  be  the  cosines  of  the  solar  and  satellite  zenith  angles, 
respectively.  Moreover,  let  the  Chs.  1  and  3  solar  reflectances  in  the  direction 
of  (/i,  A^)  for  cirrus  cloudy  conditions  be  ri(/io,M,A0)  and  , 
respectively,  which  can  be  evaluated  from  the  following  parametric  equation 
(Liou,  1980,  p.  216;  Arking  and  Childs,  1985): 


7i(M)7i(M0)’^ai 


1  -  r 


1,3. 


ai^i 


(1) 


where  r*(/io,M»^<^)  the  solar  reflectances  due  to  cloud  layer  only,  r^  are  the 

albedo  of  the  cloud  layer,  and  r^^  are  the  effective  surface  albedos  that  take 
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into  account  the  Rayleigh  scattering  effects  in  the  visible  spectral  region  and 
the  water  vapor  absorption  effects  in  the  near  IR  region.  Moreover,  7i(/i)  and 
7i(Mo)  total  directional  transraittances  in  the  direction  of  and 

(/io,0),  respectively.  These  transmittances  are  composed  of  direct  and  diffuse 
transmittances  of  the  incoming  radiation  as  follows: 

7i(M)  -  e  ^  +  tj^C/x)  , 

where  are  the  optical  depths  for  Chs .  1  and  3,  and  t^  are  given  by 

ti(M)  -  .  (3) 

where  ,ti,L4>)  are  the  downward  transmitted  diffuse  radiance  in  the  direction 

/i  for  Chs.  1  and  3,  and  Foi  are  the  respective  incident  solar  irradiances.  The 
albedo  of  clouds  can  be  obtained  from  the  following  (Liou,  1980,  p.  202): 

ri  -  2  Jq  TO^in)  /id/i,  (4) 

where  is  defined  similarly  to  Eq.  (3).  Finally,  the  second  term  in  Eq.  (1) 

represents  the  fraction  of  incoming  solar  radiation  that  is  reflected  to  the 
satellite  due  to  multiple  reflections  between  the  surface  and  the  cloud  base. 

We  develop  a  look-up  table  relating  ri(/Xo,M»A^)  based  on  Eq. 

(1)  for  different  satellite  viewing  geometries  and  cirrus  cloud  optical  and 
microphysical  properties.  In  Figure  2,  the  procedures  for  the  construction  of 
such  a  look-up  table  for  removing  the  3.7  /xra  solar  radiance  are  outlined.  We 
first  estimate  the  mean  Ch.l  effective  surface  albedo,  r^i,  using  satellite  data. 
Once  clear  pixels  are  identified  via  the  preceding  detection  scheme,  the  angular 
dependent  Ch.l  solar  reflectance,  ri(/io  i  M each  clear  pixel  is  computed 

from  the  satellite  data,  ri(/i,A^),  which  has  been  normalized  with  respect  to  the 
incident  solar  radiation  in  the  zenith  direction  as  follows: 
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Construction  of  Look-up  Table  Correlating  Chs.  1  and  3  Reflectances 


Fig  2.  A  flow  chart  for  the  construction  of  a  look-up  table  involving  the  correlation  between 
Ch.  1  and  Ch.  3  solar  reflectances. 
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(5) 


-  fi 

Subsequently,  we  construct  a  one -dimensional  histogram  of  Yi(fiQ ,  fi  for  all 
clear  pixels,  and  directly  estimate  r^^  from  the  peak  value  of  the  histogram. 
This  procedure  will  yield  reasonable  values  of  r^^,  because  the  Ch.  1  clear 
radiance  contains  the  effects  of  Rayleigh  scattering. 

The  preceding  procedure  is  not  applicable  to  evaluating  the  Ch.  3  effective 
surface  albedo,  r^s,  since  the  Ch.  3  clear  radiance  contains  both  solar  and 
thermal  infrared  components .  The  magnitude  of  the  thermal  infrared  component  is 
generally  unknown.  To  estimate  r^a,  we  use  the  following  formula  based  on  the 
radiative  transfer  principle: 

^a3  ^  ^s3 

where  is  the  total  atmospheric  transmittance  in  the  direction  of  /i,  and  rga 
is  the  Ch.  3  surface  albedo.  The  three  terms  in  Eq .  (6)  represent  downward 
atmospheric  transmission,  surface  reflection,  and  upward  atmospheric  transmission 
of  solar  radiation,  respectively.  Since  the  effects  of  Rayleigh  scattering  on 
the  transfer  of  the  3 . 7  /im  wavelength  are  relatively  weak,  multiple  reflections 
between  the  surface  and  the  Rayleigh  layer  are  neglected  in  Eq.  (6).  Based  on 
calculations  using  the  LOWTRAN  program  and  the  U.S.  Standard  Atmosphere  Profile 
(Handbook  of  Geophysics,  1985),  the  atmospheric  transmittance  between  3.6  and  4.0 
fj,m  is  on  the  order  of  0.85-0.90. 

The  exact  value  of  r53  is  also  unknown.  However,  from  laboratory 
measurements  (Sutherland,  1986)  the  3 . 7  /im  surface  emissivity  varies  between  0.8 
and  1.0  for  several  types  of  natural  soil  and  vegetation.  Thus,  the  surface 
albedo,  rg3 ,  has  a  value  between  0.0  and  0.2  based  on  Kirchhoff's  law.  Bunting 
and  Hardy  (1974)  presented  an  average  rg3  value  of  0.15  for  the  snow-free  land. 
Allen  et  al .  (1990)  reported  that  for  snow,  rs3  depends  on  the  snow  grain  size 
and  varies  between  0.005  and  0.025.  Based  on  the  preceding  values, 
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as  0.1  in  the  present  removal  scheme.  Once  the  effective  surface  albedos  for 
Chs.  1  and  3  are  determined,  we  then  proceed  to  compute  the  solar  reflectances 
of  the  two  channels  following  Eq.  (1).  For  this  purpose,  we  compute  the  single¬ 
scattering  properties  based  on  the  measured  ice  crystal  size  distributions. 
Details  of  this  step  are  described  in  the  following. 

We  have  analyzed  six  composite  ice  crystal  size  distributions  that  were 
obtained  from  the  data  presented  by  Heymsfield  and  Platt  (1984) ,  Takano  and  Liou 
(1989a),  and  the  FIRE-I-IFO  microphysic  data.  They  are  denoted  as  Cold  Ci,  - 
60°C,  Cs,  FIRE-I-IFO  Nov.  1,  FIRE-I-IFO  Nov.  2,  and  Ci  Uncinus.  We  first  define 
a  mean  effective  width  (or  size)  to  represent  ice-crystal  size  distribution  in 
the  form 


De  -  D^LnCDdL/fJ™^^  DLn(L)dL,  (7) 
J^in  J^ain 

where  D  and  L  denote  the  width  and  the  maximum  dimension  of  a  nonspherical  ice 
crystal,  respectively,  and  n(L)  is  the  size  distribution  in  terms  of  L.  The 
rationale  for  defining  to  represent  ice  crystal  size  distribution  is  that  the 
scattering  of  light  is  related  to  the  geometric  cross  section,  which  is 
proportional  to  LD.  Light  scattering  and  absorption  programs  developed  by  Takano 
and  Liou  (1989a)  for  hexagonal  ice  crystals  have  been  used  to  compute  single¬ 
scattering  properties  as  functions  of  D^.  To  compute  D^,  we  use  the  empirical 
relationship  of  L  and  D  based  on  laboratory  data  (Auer  and  Veal,  1970): 


D  -  aL*’, 


(8) 


where  a  and  b  are  empirical  coefficients.  Each  size  distribution  is  discretized 
into  five  intervals,  and  De  is  computed  by  summation  of  the  integrands  in  Eq, 

(7).  The  single-scattering  albedos  (w)  ,  phase  functions,  and  extinction 
coefficients  (^g)  ^  given  ice  crystal  size  distributions  are  computed  from 
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the  following  equations: 


w  -  n(L)dL/f^^^e  n(L)dL 

J^fiin  J^nin 


'^ax/ 


(9) 


P(e)  -  rj^ax  p(e)^gn(L)dL/fJ^^^  ;9sn(L)dL, 

J^iin  J^nin 


(10) 


n(L)dL,  (U) 

J  ^in 

where  and  are  the  scattering  and  extinction  cross  sections,  respectively, 
and  P(e)  is  the  phase  function  as  a  function  of  the  scattering  angle  6.  Table 
1  lists  the  values  of  Dg,  u),  g,  and  for  six  measured  ice  crystal  size 
distributions,  whose  g  is  the  asymmetry  , factor  associated  with  the  phase 
function.  For  0.63  /xm,  w  is  very  close  to  1  for  all  D^s ,  while  for  3.7  u> 
decreases  from  0.79  for  the  smallest  Dg  to  0.59  for  the  largest  Dg.  Moreover, 
g  for  both  wavelengths  increases  with  Dg.  The  implication  of  these  behaviors  in 
the  removal  of  the  3 . 7  /xm  solar  component  will  be  discussed  in  the  following. 

The  single -scattering  parameters  for  hexagonal  ice  columns/plates  (Takano 
and  Liou,  1989a)  are  used  in  the  adding- doubling  radiative  transfer  program 
developed  by  Takano  and  Liou  (1989b)  to  compute  r*  (/Xq  ,/x,  A(^)  ,  7i(M)  »  7i(Mo)» 
r^.  This  ice  crystal  model  is  the  best  we  have  so  far.  The  effects  of  ice 
crystal  shapes  on  the  single- scattering  properties  as  they  pertain  to  remote 
sensing  applications  is  a  subject  requiring  further  investigations.  In  the 
computations,  we  use  10  optical  depths  ranging  from  0.125  to  64.  Figure  3  shows 
an  example  of  the  relationship  between  the  reflectances  and  r^  for  the  six  ice 
crystal  size  distributions  and  the  prescribed  viewing  geometry.  Numbers  in 
parentheses  are  Dg  values.  Values  of  effective  surface  albedos  for  Chs.  1  and 
3,  x^i  and  x^^*  specified  to  be  0.12  and  0.046,  respectively.  It  is  noted 

that  radiances  of  both  channels  converge  to  the  values  of  x^i  and  x^^  for  thin 
cirrus.  Reflectance  r3  reaches  asymptotic  values  for  r^  >  0.4  due  to  saturated 
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Table  1,  Single  scattering  properties  for  six  composite  measured  size  distributions,  where  D,  is  the  mean 
effective  ice  crystal  size,  P,  is  the  extinction  coefficient,  o)  is  the  single  scattering  albedo,  and  g 
is  the  asymmetry  factor. 
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3.  Correlations  between  the  Ch.  3  solar  reflectance  (rg)  and  Ch.  1  reflectance  (r,)  for 
measured  ice  crystal  size  distributions. 


absorption  by  clouds  with  large  optical  depth.  Moreover,  this  asymptotic  value 
decreases  with  increasing  D^,  because  o)  decreases  with  increasing  and  g 
increases  with  increasing  D^.  The  preceding  analyses  have  the  following 
implications.  First,  the  removal  of  3.7  solar  component  depends  on  the  values 
of  r^i  and  r^a  for  thin  cirrus  clouds.  However,  this  dependence  becomes  less 
significant  for  r^^  >  0.4.  Thus,  it  is  important  to  specify  the  values  of  r^i  and 
r^a  as  accurate  as  possible  if  thin  cirrus  clouds  are  involved.  We  may 
improve  the  initial  guess  for  t>y  determining  x^^  the  information  of 
surface  land  type.  Second,  for  r^  >  0.4,  ra  depends  on  only.  Thus,  it  is 
possible  to  estimate  ra,  T^,  and  directly  from  Chs .  1,  3,  and  4  radiances  for 
thick  cirrus.  Finally,  for  >  100  /xm,  ra  is  independent  of  as  it  approaches 
a  minimum  value  of  --0.02  for  the  present  set  of  geometric  parameters. 

Once  the  look-up  correlation  table  is  established,  the  solar  component  in 
the  3.7  fim  radiances  is  removed  with  the  aid  of  the  0.63  fxm  reflectance  and  a 
trial  mean  effective  ice  crystal  size.  After  this  removal,  the  remaining  thermal 
infrared  radiance  of  3.7  /zm  coupled  with  the  10.9  fxm  radiance  is  used  to  perform 
retrieval  using  the  program  developed  by  Ou  et  al .  (1993)  to  obtain  the  cirrus 
cloud  temperature,  mean  effective  ice  crystal  size  (from  cloud  temperature),  and 
emissivity.  A  brief  overview  of  the  retrieval  methodology  is  given  in  the  next 
section. 

2.3  Retrieval  of  Cirrus  Cloud  Parameters 

The  retrieval  methodology  follows  the  principles  of  the  dual  thermal 
infrared  channel  technique  presented  in  Liou  et  al.  (1990)  and  Ou  et  al.  (1993). 
The  retrieval  of  cirrus  cloud  parameters  is  based  on  the  numerical  solution  of 
the  following  set  of  nonlinear  algebraic  equations  derived  from  the  theory  of 
radiative  transfer  and  parameterizations : 
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(1  ~  ^i)  Bi(T(;), 

(12) 

"  ^ai 

1  -  exp 

(-kif)  ,  i  -  3,4, 

(13) 

3 

B3(Tc) 

-  £  an[B4(Tc)]’". 

n-0 

(14) 

where  are  the  upwelling  thermal  infrared  radiances  at  the  top  of  the 
atmosphere  for  Chs .  3  and  4  over  cirrus  clouds,  are  the  corresponding  thermal 
infrared  radiances  reaching  the  cloud  base  for  the  two  channels,  are  the  cloud 
emissivities ,  is  the  equivalent  radiating  cirrus  temperature,  5^(1^)  are  the 
Planck  intensities  at  T^,  r  is  the  visible  optical  depth  (simply  referred  to  as 
optical  depth) ,  and  are  the  adjustment  factors  to  account  for  scattering  and 
wavelength  difference,  referred  to  as  the  effective  extinction  coefficients. 
Polynomial  form  is  used  to  relate  and  B^CTc)  in  Eq.  (14).  Because  both 
B3(Tc)  and  64(1^)  are  smooth  functions  of  ,  polynomial  fitting  of  these  values 
gives  very  small  errors  (<  1%)  .  The  values  of  polynomial  coefficients  are  ao  “ 
2.6327  X  10"\  ai  -  -1.063  x  10’‘\  a2  =  8.2976  x  10’^  and  ag  =  3.7311  x  10'^ 

The  solution  involves  the  effective  extinction  ratio  for  the  two  channels, 
k4/k3 ,  which  is  dependent  on  Dg.  To  investigate  the  dependence  of  the  effective 
extinction  ratio  on  Dg,  we  have  used  the  light  scattering  and  radiative  transfer 
program  developed  by  Takano  and  Liou  (1989a,  b)  to  compute  k4/k3  as  a  function 
of  Dg .  Table  2  shows  the  dependence  k4/k3  on  Dg  based  on  12  measured  size 
distributions  (Ou  et  al .  ,  1993).  The  value  of  k^/kg  increases  as  Dg  becomes 
smaller,  implying  that  more  scattering  is  associated  with  the  3.7  fim  wavelength 
than  10.9  /im.  Thus,  it  is  necessary  to  know  Dg  in  order  to  compute  k4/k3.  In 
the  dual-channel  technique  (6.5  /im  and  10.5  fim)  developed  by  Liou  et  al.  (1990), 
information  of  the  ice  crystal  size  is  not  required  because  the  dependence  of  the 
emissivity  on  the  size  is  negligible  in  these  two  wavelengths.  A  second-degree 
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polynomial  least-square  fitting  for  this  ratio  in  terms  of  l/D^  has  been 
developed: 


k4/k3  - 


2 

E 


n-o 


bn  D, 


-n 


(15) 


where  bp  *  0.722,  b^  -  55.08,  and  b2  =  174.12. 

Based  on  a  large  number  of  cirrus  microphysical  data  collected  by  optical 
probes  during  flights  over  midlatitudes,  Heymsfield  and  Platt  (1984)  have 
developed  a  functional  form  for  the  measured  ice  crystal  size  distribution  using 
a  general  power  form  as  follows: 


n(L) 


AiL®^(IWC),  L<Lo, 

A2L®2(iwC),  L>Lo, 


(16) 


where  ;  IWC  is  the  ice  water  content;  which  can  be 

parameterized  as  a  function  of  temperature  (Liou  1992)  and  A^  2  ^1,2 

temperature -dependent  coefficients  determined  from  the  measured  data.  Table  3 
lists  the  values  of  constants  in  Eq .  (16).  According  to  this  parameterization, 
Dg  is  a  function  of  temperature.  We  perform  a  least- square  polynomial  fitting 
to  relate  to  cloud  temperature  (T^.)  in  the  form 

3 

De  -  E  ‘^n(Tc  -  273)", 

n-o  ' 


where  Cq  «  326.3,  Ci  =  12.42,  C2  =  0.197,  and  C3  =  0.0012.  From  this  equation, 
Dq  increases  with  increasing  T^. ,  which  is  consistent  with  recent  aircraft 
microphysics  measurements  reported  by  Heymsfield  et  al .  (1990). 

In  order  to  solve  Eqs.  (12) -(14),  the  upwelling  radiances  reaching  the 
cloud  base,  must  be  given.  They  are  determined  by  using  the  radiances  of 
clear  pixels,  because  very  little  water  vapor  is  present  above  cirrus  clouds. 
It  follows  that  if  we  can  identify  the  clear  pixels,  may  be  determined  by  a 
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Table  2  Values  of  the  effective  extinction  ratio  k^/ks  as  a  function  of  the 
mean  effective  size  D^.  Values  of  are  computed  from  a  number  of  measured  ice- 
crystal  size  distributions. 


k*/k3 

23.9 

2.725 

30.4 

2.364 

41.5 

1.792 

47.5 

1.908 

57.9 

1.705 

64.1 

1.585 

75.1 

1.347 

93.0 

1.231 

104.1 

1.242 

106.7 

1,224 

110.4 

1.212 

123.6 

1.182 

statistical  method.  We  may  select  a  scene  and  use  the  data  within  this  scene  to 
construct  a  two-dimensional  histogram  in  the  domain  of  R3  and  .  The  radiances 
corresponding  to  the  peak  of  the  frequency  distribution  are  assigned  as  the  mean 
clear  radiances.  The  area  of  the  scene  should  be  large  enough  to  contain  a 
statistically  significant  amount  of  pixels.  At  the  same  time,  this  area  should 
also  be  sufficiently  small  to  guarantee  the  homogeneity  of  the  surface 
temperature  and  water  vapor  distributions  within  the  scene.  Ou  et  al.  (1993) 
used  the  FIRE-I-IFO  data  and  two  sizes  (1*"  x  1°)  and  (0.5°  x  0.5°)  to  investigate 
the  effects  of  varying  the  area  size  on  the  retrieval  and  found  that  the 
variability  of  the  clear  radiances  decreases  with  decreasing  size  of  the  area. 
However,  the  mean  clear  radiances  remain  unchanged.  Both  Minnis  et  al.  (1993b) 
and  Chou  (1991)  used  the  ISCCP  grid  (2.5°  x  2.5°)  for  the  determination  of  clear 
radiances . 

2.4  Removal -Retrieval  Scheme 

Both  the  removal  and  retrieval  programs  are  dependent  on  the  ice  crystal 
size  distribution  represented  by  the  mean  effective  size,  Dg.  The  mean  effective 
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Table  3.  Values  of  A^,  A2,  Bi,  B2,  Lq,  and  IWC  in  Eq.  (16)  as  functions  of 
temperature . 


Temperature 

(°C) 

-20  - 

-25 

-25  - 

-30 

-30  - 

-35 

-35  - 

-40 

-40  - 

-45 

-45  - 

-50 

-50  - 

-55 

-55  - 

-60 

6.837x10®  3.304x10^2 


7.328x10®  3.068x10^* 


1.954x10®  9.063x10^2 


7.540x10®  1.329x10“ 


2.163x1012 


1.117x10“ 


1.778x10“ 


2.800x10“ 


IWC 

(g/m®) 


0.0270 


0.0250 


0.0175 


0.0126 


0.0034 


0.0025 


0.0018 


0.0009 


sizes  are  computed  on  a  pixel-by-pixel  basis.  They  are  determined  from  the 
retrieved  cloud  temperature  based  on  the  relationship  between  the  mean  effective 
size  and  cloud  temperature  (Eq.  17)  and  not  selected  out  of  the  six  measured  size 
distributions  presented  in  Figure  3.  A  value  of  is  prescribed  initially. 
This  value  is  usually  the  largest  possible  value  of  D^,  which  corresponds  to  the 
minimum  amount  of  solar  radiances  to  be  removed.  Using  satellite-derived  r^i, 
and  prescribed  rs3,  a  look-up  table  for  r3  versus  r^  is  constructed  according  to 
the  procedures  described  in  Section  2.2.  The  solar  reflectance  r3  is  then 
obtained  from  the  measured  0.63  /im  reflectance  and  the  look-up  table.  Let  the 
reflected  radiance  of  Ch.  3  be  Rssoi-  Thus, 

R3so1  -  /‘0F03r3. 


where  F03  is  the  incident  solar  irradiance,  which  is  obtained  from  Thekaekara 
(1974),  and  r3  is  the  solar  reflectance  for  the  3.7  fim  channel.  Let  the  thermal 
infrared  component  of  Ch.  3  be  initial  value  for  R^ir  obtained  from 
the  following 
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(19) 


^3ir 


(0)  _ 


R3 


Ro 

^3sol 


Both  R3ir^°^  and  R4,  [Eq.  (12)]  are  used  in  the  cirrus  retrieval  program. 
The  same  as  in  the  removal  scheme  is  utilized  to  specify  k4/k3  according  to 

Eq.  (15).  A  new  cirrus  cloud  temperature  is  then  obtained.  From  we 

may  compute  an  intermediate  value  for  using  Eq.  (17).  Because  we  have 

removed  the  minimum  amount  of  the  3.7  /im  solar  component,  the  thermal  infrared 
component  Rair^^^  will  be  overestimated  based  on  Eq.  (19).  From  Eqs .  (12)  and 
(17),  and  will  be  underestimated.  Thus,  the  difference  ~ 

would  be  substantial.  If  we  apply  to  the  removal -retrieval  scheme, 

numerical  instability  may  occur.  To  avoid  this,  a  new  suppressed  value  of 
is  computed  as  follows: 


(20) 


The  suppression  factor  a  is  on  the  order  of  0.005.  The  use  of  such  a  small  value 
for  a  maintains  the  stability  in  numerical  iterations.  Using  the  new  value  for 
we  repeat  the  above  removal -retrieval  procedures  to  obtain  The 

iteration  process  is  carried  out  until  satisfies  the  following  convergence 

criterion: 


|^(n)  .  p(n-l)' 


<  1.0x10"*^. 


(21) 
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3. 


SENSITIVITY  STUDIES  AND  ERROR  ANALYSES 


Sensitivity  studies  have  been  performed  using  the  midlatitude  summer 
atmospheric  profile.  First,  we  must  generate  simulated  thermal  infrared  and 
solar  clear  and  cloudy  radiances.  Spectral  transmittances  for  the  3.7  /im  and 
10.9  /im  wavelengths  have  been  computed  by  using  the  program  developed  by  Weinreb 
and  Hill  (1980),  in  which  the  instrument  filter  function  is  accounted  for.  The 
simulated  upwelling  IR  radiances  for  clear  sky  are  computed  using  the  solution 
of  the  basic  radiative  transfer  equation  in  the  forms 

Ra3ir  -  B3(Ts)T3(ps)  H- JJ®B3(T)£I|Mdp,  (22a) 

Ra4  -  B4(Ts)T4(Ps)  +  64(1)  dp .  (22b) 

where  Pg  and  T^  are  the  surface  pressure  and  temperature,  respectively,  and  T^Cp) 
are  the  atmospheric  transmittances  at  pressure  level  p.  In  the  simulations, 
cloud-base  heights  are  prescribed  at  7 ,  9,  and  11  km.  The  cloud  temperatures  and 
mean  ice  crystal  sizes  can  be  determined  from  the  temperature  profile.  For  each 
cloud-base  height,  calculations  are  carried  out  for  10  optical  depths  (0.125, 
0.25,  0.5,  1.0,  2.0,  4.0,  8.0,  16.0,  32.0,  64.0).  Emissivity  values  are  obtained 
from  Eq.  (13)  .  The  upwelling  thermal  infrared  cloudy  radiances  are  then  computed 
from  Eq.  (12) . 

To  compute  the  Ch.  3  solar  radiance  for  clear  and  cloudy  skies,  we  first 
set  the  effective  surface  albedos  for  Chs .  1  and  3  to  be  0.12  and  0.1, 
respectively.  The  geometric  parameters  used  are:  solar  zenith  angle,  6q  «  71®; 
satellite  zenith  angle,  6  ■=  40®;  and  relative  azimuthal  angle,  =  146®.  This 
set  is  based  on  the  geometric  parameters  involving  the  FIRE-I-IFO  data.  The 
clear  solar  radiance  for  Ch.  3  is  computed  as  follows: 
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Figure  4(a)  and  4(b)  show  the  distributions  of  and  R3  (“  Rssoi  ^air) 
versus  R^,,  respectively,  for  the  three  cloud  heights  and  10  optical  depths.  In 
each  plot,  the  Planck  intensity  curve  is  included  as  a  reference.  The  curves  of 
constant  cloud  height  are  all  above  the  Planck  intensity  curve.  The  uppermost 
point  where  all  the  curves  converge  represents  clear  radiances.  The  lowest  point 
for  each  curve  corresponds  to  the  highest  optical  depth.  The  curvature  of  the 
curves  in  Figure  4(b)  is  due  to  the  variation  in  Rasoi*  comparing  Figures  4(a) 
and  4(b),  we  note  that  R3  is  larger  than  ^Y  ^t)out  0.05-0.2  mW/mVsr/cm*’^  due 
to  additional  solar  reflections  for  the  present  sun-satellite  geometry.  This 
difference  will  increase  as  the  sun  becomes  closer  to  the  zenith.  Moreover,  in 
Figure  4(a),  when  the  emissivity  approaches  1  for  the  largest  optical  depth,  both 
R3ir  and  R^  approach  and  54,(1^),  respectively.  However,  in  Figure  4(b),  the 
total  Ch.  3  radiance  R3  >  ,  The  difference  between  R3  and  B3(Tc)  is  mainly 
due  to  the  size  effect  on  the  solar  component  as  well  as  the  atmospheric 
absorption  and  surface  reflection. 

Figure  5  illustrates  the  ratio  of  R3S0I  ^3ir  ^  function  of  r  and  cloud- 
base  height.  Variation  in  the  relative  magnitude  of  R3S01  ^3ir  seen. 
For  the  present  sun- satellite  geometry,  this  ratio  increases  with  r.  For  thin 
cirrus  (r  <  1),  Rg^^  is  10  times  larger  than  R3S0I  is  almost  independent  of 
cloud  height.  However,  for  thick  cirrus,  RssoiAsir  approach  asymptotic  values, 
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(a)  Midlatitude  summer 

•-  Cloud  height  =  7  km 

Cloud  height  =  9  km 

Cloud  height  =  11  km 

Planck  intensity  curve 
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(b)  Midlatitude  summer 
•-  Cloud  heights  7  km 
0-  Cloud  heights  9  km 
■-  Cloud  height  =  1 1  km 
—  Planck  intensity  curve 


0^  =  71°,  e  =  40^A<|)  =146= 

'■ai  =  0  ‘'2,  ra3  =  0.1 
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(a)  Distribution  of  the  simulated  Ch.  3  thermal  infrared  (Rgj^)  and  Ch.  4  (R^)  radiances  for  three 
cloud  heights.  Also  depicted  in  the  figure  is  the  Planck  intensity  curve,  (b)  Distribution  of  the 
simulated  Ch.  3  solar  plus  thermal  infrared  (Rg)  and  Ch.  4  (R^)  radiances  for  the  same  three 
cloud  heights. 
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because  R^ir  approaches  the  Planck  intensity  value  at  the  cirrus  cloud 
temperature.  Rssoi  also  approaches  its  asymptotic  value.  Moreover,  the  radiance 
ratio  increases  with  the  cloud-base  height,  because  of  the  coupling  of  lower 
cloud  temperatures  and  smaller  ice  crystals.  The  former  leads  to  smaller  R^h 
values,  while  the  latter  produces  stronger  solar  reflection.  For  the  present 
sun-satellite  geometry,  Figure  5  shows  that  the  Ch.  3  solar  component  is  at  least 
10%  of  the  total  radiance  which  can  not  be  ignored. 

In  the  following,  we  investigate  the  effects  of  uncertainties  in  the  input 
parameters  on  the  accuracy  of  the  removed  solar  component.  Errors  in  the  3.7  fim 
solar  component  can  be  due  to  uncertainties  in  the  mean  effective  size  (D^)  ,  as 
well  as  the  Chs .  1  and  3  effective  surface  albedos  (r^i  and  . 

First,  the  effects  of  uncertainties  in  on  the  accuracy  of  the  Ch.  3 
solar  radiance  (R3S01)  studied.  From  the  three  prescribed  cloud-base  heights 
of  7,  9,  and  11  km,  mean  values  of  136,  81,  and  52  /xm  are  obtained.  Next, 
^3soi  total  Ch.  3  radiance  (R3)  are  computed  for  each  value,  as  well 
as  for  a  number  of  optical  depths.  We  assume  that  r^i  and  exact 
values  and  that  the  value  of  used  in  the  scheme  has  a  range  of  uncertainty, 
ADq  -  ±5  /xm.  We  then  proceed  to  quantify  R3S01  by  using  the  look-up  table,  as 
displayed  in  Figure  3.  Errors  in  R3S0I  will  result  due  to  the  uncertainty  in  D^. 
Figure  6(a)  shows  the  maximum  percentage  errors  for  AR^^^i/R^^qi  as  functions  of 
optical  depth.  The  relative  error  AR3soi/R3soi^  increases  with  the  optical  depth. 
For  small  optical  depths,  AR^^^i  is  small  because  R3S01  is  dominated  by  surface 
reflection.  For  thicker  clouds,  the  reflected  radiances  are  dominated  by  cloud 
reflection  and  the  effects  of  ADg  become  important.  Since  R3S0I  approaches  an 
as3rmptotic  value  for  R^  >  0.4  (Figure  3)  and  since  the  range  of  AD^  is  fixed  at 
±5  /xm,  AR3S01/R3S0I  also  approaches  an  asymptotic  value.  The  error  in  AR3S01/R3S0I 
for  the  7 -km  cloud  case  is  smaller  than  the  errors  for  the  9-  and  11 -km  cloud 
cases,  because  the  former  has  a  larger  than  100  /xm.  Note  that  R3S0I  is 
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(a)  Midlatitude  summer 
0Q  =  7r,  0  =  40°,  A<l)  =  146® 

fai  =  0-12.  r,3  =  0.1 


AD-  =  ±5  pm 


Cloud  height  =  7  km 
Cloud  height  =  9  km 


Hi-  Cloud  height  =  1 1  km 


Ootical  Depth 


Optical  Depth 

Fig  6.  Maximum  percentage  errors  in  the  Ch.  3  solar  radiances  (R3gQ|)  produced  by  the  errors  in  (a) 
the  mean  effective  size.  AD^,  of  ±  5  pm,  (b)  the  Ch.  1  mean  effective  surface  albedo,  Ar^^ ,  of 
±0.02,  (c)  the  Ch.  3  mean  effective  surface  albedo,  Ar  o,  of  ±  0.05. 


nearly  constant  for  >  100  fim  (Figure  3).  Overall,  the  maximum  error  of 
^Rasoi/Rssoi  is  less  than  12%. 

Next,  we  investigate  the  effects  of  uncertainties  in  r^^  on  the  accuracy 
of  Rssoi-  We  assume  that  and  r^s  are  both  known  without  errors,  but  r^^  has  an 
uncertainty  of  ±0.02.  Errors  in  R^soi  >^ill  be  produced  due  to  the  uncertainty  in 
r^i.  Figure  6(b)  shows  the  maximum  percentage  errors  in  Rasoi  ^s  functions  of 
optical  depth.  In  contrast  to  Figure  6(a),  AR3S01/R3S01  decreases  as  the  optical 
depth  increases.  This  decrease  occurs  because  the  reflected  radiances  are 
essentially  controlled  by  cloud  reflection  for  larger  optical  depths,  so  that  the 
changes  in  the  effective  surface  albedo  have  little  effect.  For  smaller  optical 
depths,  AR3S01/R3S03L  is  larger  for  the  7 -km  cloud  than  for  the  9-  and  11-km  clouds, 
because  R3S0I  is  much  smaller  for  larger  D^s  associated  with  the  7 -km  cloud. 
Overall,  for  optical  depths  larger  than  1,  the  maximum  AR3so]^/R3so]^  is  less  than 
12%. 

We  also  investigate  the  effects  of  uncertainties  in  on  the  accuracy  of 
^Ssoi*  is  assumed  that  and  r^i  values  are  exact.  However,  r^a  has  an 
uncertainty  of  ±0.05,  which  will  result  in  errors  in  Rasoi-  Figure  6(c)  shows  the 
maximum  percentage  errors  in  R3S0I  ^s  functions  of  optical  depth.  As  in  Figure 
6(b),  AR3S01/R3S01  decreases  as  the  optical  depth  increases.  In  addition,  for 
smaller  optical  depths,  is  smaller  for  the  11-km  cloud,  because  the 
value  of  Rssoi  is  larger  for  smaller  associated  with  higher  clouds.  For 
optical  depths  larger  than  1,  the  maximum  AR3S01/R3S0I  i®  less  than  15%. 

Errors  in  the  retrieved  cirrus  parameters  can  be  due  to  instrument  noises 
and  the  errors  in  r^x  and  r^s.  To  perform  error  analyses,  the  probability  of  the 
error  produced  by  a  particular  parameter  x»  Ax,  is  assumed  to  follow  the  Gaussian 
distribution  in  the  form: 
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P(Ax)  - 


1 


(26) 


_ _  exp[-(Ax^/2a2) ] , 

where  a  is  the  standard  deviation.  The  Gaussian  distribution  differs  from  random 
distribution  in  that  the  latter  would  give  a  uniform  probability  of  error  within 
the  prescribed  error  bounds,  while  the  former  would  give  a  decreasing  probability 
of  error  as  the  error  increases.  The  Gaussian  distribution  appears  to  be  more 
representative  of  a  realistic  error  distribution  than  the  random  distribution. 
We  have  used  an  established  software  as  the  Gaussian  distribution  generator  to 
simulate  various  error  distributions  with  the  following  specified  standard 
deviations:  for  the  Ch.  3  noise,  a  ^  0.4  K;  for  the  Ch.  4  noise,  a  —  0.12  K;  for 
the  Ch.  1  effective  surface  albedo,  a  =  0.02;  and  for  the  Ch.  3  effective  surface 
albedo,  a  *  0.05.  The  a  values  for  Ch.  3  and  Ch.  4  noises  are  prescribed  based 
on  the  values  described  by  Dudhia  (1989).  The  a  value  for  the  Ch,  1  surface 
albedo  is  based  on  the  realistic  histogram  for  Ch.  1  reflectances.  The  a  value 
for  the  Ch,  3  mean  effective  surface  albedo  is  based  on  the  range  of  laboratory 
measurements  (Sutherland,  1986)  for  3.7  fim  surface  reflectances. 

To  carry  out  the  error  analysis,  we  generate  three  thousand  sets  of  errors 
in  both  instrument  noises  and  the  prescribed  parameters.  Figures  7  (a) -(d)  show, 
respectively,  the  root-mean- square  (rras)  errors  for  the  Ch.  3  solar  component, 
the  cirrus  cloud  temperature,  the  mean  effective  size,  and  the  optical  depth,  as 
functions  of  the  prescribed  optical  depth  for  three  cloud-base  heights.  In  all 
cases,  rms  errors  decrease  with  increasing  optical  depths  because  the  effects  of 
the  uncertainty  in  effective  surface  albedo  on  cloudy  radiances  become  less  as 
the  optical  depth  increases.  Moreover,  for  thin  cirrus  clouds,  slight 
perturbations  in  both  radiances  due  to  instrument  noises  would  lead  to  large 
errors  in  the  retrieved  parameters.  For  optical  depths  larger  than  8,  the 
percentage  error  in  the  optical  depth  increases  [Figure  7(d)]  because  the 
emissivity  approaches  1  and  a  slight  error  in  it  would  lead  to  a  large  error  in 
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in  Ch.  3  surface  albedo  (a  =  0.05)  as  functions  of  the  optical  depth. 


the  optical  depth.  For  optical  depths  greater  than  0.25,  the  rras  errors  in 
Tj.,  Dq,  and  r,  are  less  than  0.0045  raW/mVsr/cm‘^  (or  5%),  2.6  K,  15  /xm,  and  6%, 
respectively,  indicating  that  the  removal -  retrieval  scheme  is  numerically  stable 
and  the  retrieval  accuracy  is  acceptable.  The  rms  errors  in  are  of  the  same 
order  of  magnitude  as  the  rms  differences  between  the  retrieved  and  lidar  derived 
values  (Table  3,  Minnis  et  al.,  1993b).  Moreover,  the  rms  errors  in  are  less 
than  the  uncertainty  in  the  determined  from  the  aircraft  measurements 

(Heymsfield  and  Platt,  1984). 

The  retrieval  accuracy  depends  on  the  cloud  height  and  varies  for  the 
retrieved  cloud  parameters.  Errors  in  Rasoi  increase  with  decreasing  cloud 
height.  Errors  in  T^  are  virtually  independent  of  cloud  height  for  optical 
depths  greater  than  2.  However,  for  optical  depths  smaller  than  2,  errors 
increase  with  decreasing  cloud  height.  This  behavior  is  related  to  errors  in  Dg, 
which  increase  with  decreasing  cloud  height,  as  shown  in  Figure  7(c).  In  the 
retrieval  scheme,  the  parameters  D^,  k^/k3,  and  T^  and  the  associated  errors  are 
interrelated  through  Eqs .  (15)  and  (17).  Moreover,  when  the  cloud  is  lower, 

cloud  temperature  is  higher,  while  ice  crystal  size  is  larger  and  k4/k3 
approaches  its  asymptotic  values  in  the  case  of  lower  clouds.  A  small  error  in 
k^/ka  would  lead  to  large  errors  in  both  Dg  and  T^ .  Finally,  in  Figure  7(d), 
errors  in  the  optical  depth  are  shown  to  increase  with  decreasing  cloud  height. 
This  implies  that  errors  in  the  related  cloud  emissivity  increase  with  cloud 
temperature . 

Effects  of  other  uncertainties  on  the  accuracy  of  the  removal -retrieval 
scheme  are  discussed  in  the  following.  Based  on  Figure  3,  errors  in  the  Ch.  1 
radiance  do  not  affect  significantly  on  the  accuracy  of  the  removed  solar 
component  for  optical  depths  larger  than  one,  because  the  Ch.  3  radiance  is  close 
to  the  asymptotic  value.  For  optical  depths  smaller  than  1,  errors  of  the 
removed  solar  component  depend  on  the  effective  surface  albedo  for  both  channels 
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and  the  ice  crystal  sizes.  The  inhomogeneity  of  cirrus  clouds  would  have  an 
important  effect  on  their  radiative  properties.  However,  it  is  not  feasible  to 
retrieve  the  detailed  structure  of  clouds  utilizing  the  infrared  methodology. 
The  retrieval  results  would  represent  the  averaged  cloud  characteristics  over  the 
AVHRR  pixels.  It  appears  to  us  that  the  effects  of  ice  crystal  shape  and  cloud 
inhomogeneity  on  the  radiative  properties  of  clouds  will  require  considerable 
research  before  their  impacts  on  satellite  remote  sensing  can  be  understood. 
Finally,  regarding  the  effects  of  errors  in  the  polynomial  fitting  of  k^/ks  vs. 
De,  we  have  shown  that  the  root -mean- square  difference  between  the  computed  and 
least-squares  fitted  values  of  k^/ks  is  0.045,  revealing  that  the  polynomial 
fitting  is  an  excellent  parameterization  to  relate  k^/ks  and  the  mean  effective 
size  (Ou  et  al.,  1993).  Thus,  it  is  not  necessary  to  study  on  the  sensitivity 
of  errors  in  the  parameterization  of  k^/k3  on  the  retrieval. 
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4. 


APPLICATION  OF  SATELLITE  DATA  ASSOCIATED  WITH  FIRE-I-IFO 


For  the  purpose  of  testing  the  removal -retrieval  scheme,  we  have  acquired 
the  AVHRR  Global  Area  Coverage  (GAC)  data  collected  by  NOAA  9  at  about  2100  UTC, 
28  October  1986.  This  data  set,  including  the  FIRE-I-IFO  region  (42°-47®  N;  87°- 
92®  W) ,  covers  the  midwestern  United  States.  It  contains  the  scaled  values  of 
latitude,  longitude,  solar  zenith  angle,  the  zenith-normalized  Chs .  1  and  2 
radiances,  and  the  brightness  temperatures  of  Chs.  3,  4,  and  5  for  each  GAC 
pixel.  There  are  a  total  of  2.6  x  10^  pixels.  A  program  has  been  developed  to 
unpack  the  data  set,  using  the  calibration  procedures  documented  by  Kidwell 
(1991)  to  obtain  the  true  data  values.  We  have  also  developed  a  program  to 
compute  the  satellite  zenith  angle  and  the  relative  azimuthal  angle  based  on  the 
viewing  geometry  of  the  polar -orbiting  satellite.  The  thermal  infrared  channels 
(3,  4,  and  5)  are  constantly  calibrated  by  blackbody  calibrators  on  board  the 
satellite.  However,  there  are  no  in-flight  calibration  facilities  to  monitor  the 
long-term  stability  of  visible  channels  (1  and  2).  Staylor  (1990)  has  estimated 
that  the  rate  of  degradation  for  these  visible  channels  is  about  6%  per  year. 
We  have  accounted  for  these  effects  in  the  calibration  procedures. 

Using  the  foregoing  data,  we  have  performed  retrievals  of  cirrus  cloud 
parameters  over  three  areas  with  different  scales,  as  shown  in  Figure  8.  We 
first  apply  the  removal -retrieval  scheme  to  a  1®  x  1®  area  over  western  Wisconsin 
(the  innermost  frame  in  Figure  8) .  We  then  expand  the  domain  of  retrievals  to 
a  5'*  X  5®  area  associated  with  the  FIRE-I-IFO  region  (the  intermediate  frame  in 
Figure  8).  These  retrievals  are  designed  for  developing  and  testing  an  automated 
removal -retrieval  scheme,  for  verifying  the  retrieval  results  against  local  lidar 
measurements,  and  for  comparing  the  retrieval  results  with  those  determined  from 
different  satellite  data  sets  and  retrieval  methodologies.  The  third  retrieval 
domain  is  the  20®  x  30®  area  over  the  midwestern  United  States.  These  retrievals 
enable  us  to  test  the  removal -retrieval  scheme  over  a  large  area  with  different 
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30°  region  covering  the  midwestern  United  States. 


satellite  viewing  geometries  and  surface  characteristics.  We  also  compare  our 
retrieval  results  with  the  AVHRR  visible  cloud  picture  for  the  same  time  period 
and  area. 

4.1  Retrievals  over  a  1®  x  1®  Area 

We  have  tested  the  removal -retrieval  scheme  over  a  small  area  (1®  x  1®) 
west  of  Fort  McCoy.  Results  of  this  test  are  presented  in  the  following.  The 
selected  zone  extends  from  43.5®  to  44.5®N  and  from  91®  to  92®W.  Within  this 
area,  there  are  392  GAG  pixels.  Figures  9 (a) -(c)  show  the  halftone  display  of 
the  Ch.  1  reflectances  and  the  Ch.  3  and  Ch.  4  brightness  temperatures.  The 
darker  areas  in  each  frame  correspond  to  lower  reflectances  or  higher  brightness 
temperatures,  and  are  associated  with  either  clear  or  very  thin  cirrus  clouds. 
The  whiter  areas  correspond  to  higher  reflectances  or  lower  brightness 
temperatures,  and  are  associated  with  thicker  cirrus  clouds.  The  area  is  mostly 
covered  with  clouds  during  the  day.  Small  clear  areas  are  evident  over  the 

eastern  and  southwestern  regions.  High  reflectances  and  low  brightness 

temperatures  preside  over  the  south- central  region.  Brightness  temperatures  are 
distributed  more  evenly  over  the  northern  half  of  the  area  than  those  over  the 
southern  half.  This  suggests  that  at  the  time  of  satellite  overpass,  the  air 
mass  over  the  southern  half  was  more  dynamically  active  than  the  northern  half. 
The  clear  detection  scheme  is  applied  to  the  satellite  data  covering  the 
aforementioned  area.  Figure  10  shows  the  distribution  of  clear  (•)  and  cloudy 
(0)  pixels.  There  are  46  detected  clear  pixels,  which  are  scattered  over  the 
eastern  and  southwestern  regions.  The  clear  areas  shown  in  Figure  9(a)  and  9(b) 
are  correctly  identified  in  Figure  10. 

To  remove  the  3 . 7  /xm  solar  component,  we  must  first  determine  r^i  and  ra3. 
The  Ch.  1  effective  surface  albedo  (-0.12)  is  taken  to  be  the  mode  value  of  r^ 
for  all  the  identified  clear  pixels.  The  determination  of  is  somewhat  more 
involved.  As  described  earlier,  the  surface  albedo  r53  is  fixed  at  0.1. 
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Fig  9.  Half-tone  display  of  (a)  the  Ch.  1  reflectance,  (b)  the  Ch.  3  brightness  temperature,  and 
(c)  the  Ch.  4  brightness  temperature,  for  the  1®  x  1°  area  ( 43.5®  -  44.5®N,  91®  -  92®W) 
over  Wisconsin  at  2100  UTC,  28  October  1986. 
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Fig  1 0.  Distribution  of  the  dear  and  doudy  pixels  that  are  identified 
over  the  1®  X  1°  area  west  of  Fort  McCoy. 


However,  to  determine  r^s  from  Eq.  (6),  the  value  of  is  required.  Two  sets 
of  sounding  data  over  the  region  at  2100  UTC  are  used  in  the  calculations.  One 
set  is  the  Green  Bay  sounding  presented  by  Starr  and  Wylie  (1990),  while  the 
other  is  derived  from  the  ISCCP-Cl  sounding  profiles  (Rossow  and  Schiffer,  1991). 
The  ISCCP  sounding  profiles  are  the  averages  of  satellite-derived  sounding 
profiles  over  a  2.5®  x  2.5®  area  at  nine  vertical  levels.  Generally,  the 
temperatures  of  these  two  profiles  differ  by  less  than  10  K,  while  the  humidities 
differ  by  less  than  2g/kg  above  700mb.  The  transmittance  profiles  for  the  3.7 
fim  wavelength  are  computed  using  the  two  sets  of  soundings.  The  total 
transmittance  in  both  cases  is  about  0.99  because  of  the  low  humidity  content, 
implying  that  the  simulated  Ch.  3  transmittance  is  relatively  insensitive  to  the 
uncertainty  in  sounding  profiles . 

We  apply  the  removal -retrieval  scheme  to  each  cloudy  pixel  that  has  been 
identified.  As  investigated  in  Ou  et  al .  (1993),  the  cirrus  cloudy  condition  can 
be  found  as  the  brightness  temperature  difference  between  Ch.  3  (nighttime)  and 
Ch.  4  larger  than  2K.  The  Ch.  3  daytime  data,  which  is  contaminated  by  the  solar 
reflection,  cannot  be  used  to  identify  the  presence  of  cirrus.  However,  cirrus 
pixels  can  be  detected  from  the  retrieved  cloud  temperature.  Figures  11 (a) -(c) 
show  the  halftone  displays  of  the  retrieved  T^,  D^,  and  r  over  the  1®  x  1®  area, 
respectively.  The  distribution  of  T^  over  the  northwestern  region  is  relatively 
homogeneous  with  values  between  230  and  240  K.  Sharp- transition  regions  are 
evident  near  cold  cirrus  clouds  (T^  <  225  K)  .  These  features  result  from  the 
presence  of  cold  cirrus  within  a  clear  region.  The  cold  cirrus  clouds  appear  to 
be  associated  with  cirrus  generation  regions  (Heymsfield,  1975) .  These  clouds 
are  relatively  high  and  thin  and  contain  small  ice  crystals  [see  Figure  11(b)  and 
(c)]. 

From  the  distribution  of  cirrus  cloud  temperature,  we  may  also 
qualitatively  infer  the  dynamic  activity  of  the  air  mass.  To  the  north  and  east 
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1 .  Half-tone  display  of  the  (a)  retrieved  cbud  temperature,  (b)  the  retrieved 
mean  effective  size,  and  (c)  the  retrieved  optical  depth,  over  the  1°  x1°  area 
west  of  Fort  Mcloy  at  2100  UTC,  28  October  1 986. 
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of  the  sharp- transition  regions,  the  distribution  of  cirrus  cloud  temperatures 
is  less  variant.  A  transition  is  seen  from  the  cirrus  growing  regions,  located 
at  the  southwestern  side  of  the  cloudy  area,  to  the  cirrus  dissipating  regions, 
located  at  the  northeastern  side.  Considering  that  the  prevalent  wind  is  from 
the  west-south-west  direction  (Starr  and  Wylie,  1990),  this  transition  pattern 
may  be  explained  by  the  theory  of  cirrus  formation-dissipation  proposed  by 
Heymsfield  (1975)  .  The  stable  regions  over  the  northern  part  show  that  the 
cirrus  clouds  there  are  dissipating.  We  conclude  that  the  cirrus  growing  regions 
on  the  southwestern  side  are  dynamically  active  with  vertical  motions  enhanced 
by  latent  heat  release,  whereas  the  cirrus  dissipating  regions  on  the 
northeastern  side  are  dynamically  stable.  It  is  likely  that  within  a  small  area 
on  the  order  of  100  x  100  km^,  newly  formed  cirrus  coexists  with  mature  and 
dissipating  cirrus. 

In  Figure  11(b),  the  distribution  of  corresponds  to  the  distribution  of 
Tq.  Sharp- transition  regions  occur  over  the  southern  and  eastern  parts  of  the 
area.  Over  the  northern  part,  the  distribution  of  De  is  relatively  uniform,  the 
distribution  of  r,  shown  in  Figure  11(c),  is  also  relatively  uniform  (r  0  - 
2.5),  except  near  the  south-central  part,  where  r  reaches  a  value  greater  than 
5.  It  is  noted  that  the  distribution  of  r  does  not  necessarily  match  the 
distribution  of  T^.  Most  of  the  cold  cirrus  and  warm  cirrus  regions  show 
relatively  small  values  of  r  (0-1.5).  Thick  cirrus  clouds  occurring  over  the 
south-central  part  represent  a  region  with  high  ice  water  path. 

4.2  Retrievals  over  a  5°  x  5°  Area 

We  have  also  applied  the  removal -  retrieval  scheme  to  the  AVHRR  daytime 
data,  collected  on  28  October  1986  at  2100  UTC  over  Wisconsin  (42®  to  47®N;  87®W 
to  92 ®W)  during  the  FIRE-I-IFO.  The  brightness  temperature  of  Ch.  4  over  this 
area  is  shown  in  Figure  12.  The  majority  of  the  cloudy  areas  is  located  near  the 
northwest  corner,  at  about  44®N  and  91®W  and  45®N  and  89®W.  These  areas  are  all 
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characterized  by  brightness  temperatures  lower  than  260  K.  Clear  areas  are 
located  mainly  at  45®N  and  90®W  and  near  the  northeast  and  southeast  corners. 

To  reduce  the  effects  of  surface  inhomogeneity  on  the  mean  clear  radiances, 
the  whole  5®  x  5®  area  is  divided  into  25  1®  x  1®  boxes  to  which  the  clear/cloud 
detection  scheme  is  applied.  Using  the  brightness  temperatures  of  clear  pixels 
that  are  identified,  we  compute  a  pair  of  Ch.  3  and  Ch.  4  thermal  infrared  mean 
clear  radiances  for  each  box.  Because  of  the  required  criteria  for  detecting 
clear  pixels,  some  boxes  contain  no  identifiable  clear  pixels.  In  these  cases, 
the  mean  clear  radiances  of  a  neighboring  box  are  used.  The  removal -retrieval 
scheme  is  then  applied  to  each  box. 

Figures  13(a) -(d)  show  the  geographic  distribution  of  the  retrieved  cloud 
temperature,  mean  effective  ice  crystal  size,  optical  depth,  and  cloud  height 
over  the  5®  x  5®  area,  respectively.  In  Figure  13(a),  the  retrieved  cloud 
temperatures  range  between  215  and  240  K.  The  distribution  of  the  retrieved 
cloud  temperature  is  different  from  that  of  the  Ch .  4  brightness  temperatures. 
For  example,  near  the  northwest  corner,  where  the  Ch.  4  brightness  temperatures 
are  about  240  K,  the  minimum  cloud  temperature  is  about  225  K.  Around  45®N  and 
89 ®W,  the  Ch.  4  brightness  temperatures  are  also  about  240  K,  while  the  cloud 
temperatures  are  about  235  K.  However,  one  degree  longitude  east  of  this  cloudy 
area,  the  Ch.  4  brightness  temperatures  are  about  255  K,  while  the  minimum  cloud 
temperature  is  about  220  K.  It  is  possible  that  cirrus  clouds  are  on  top  of  low 
clouds  near  45®N  and  89®W.  However,  a  single  layer  of  high,  thin  cirrus  clouds 
is  present  near  45®N  and  88®W.  Over  the  region  where  multi-layered  clouds  are 
present,  retrievals  had  been  performed  as  though  there  were  only  a  single  layer 
of  cloud.  This  is  because  reliable  methods  are  yet  to  be  developed  for  detecting 
multi-layer  cloudy  pixels  using  AVHRR  data.  In  Figure  13(b),  the  distribution 
of  the  retrieved  mean  effective  ice  crystal  sizes,  which  range  between  40  and  100 
/ira,  are  in  agreement  with  that  of  the  cloud  temperature  shown  in  Figure  13(a). 
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Ch.  4  Brightness  Temperature  (K) 


Fig  12.  Half  -  tone  display  of  the  Ch.  4  brightness  temperature  for  the  5°  x  5®  area  over  Wisconsin 
at  2100  UTC,  28  October  1986.  The  symbols  "+"  identify  the  location  of  three  FIRE-I-IFO 
stations  (see  text). 
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Longitude  (“W)  Longitude  (“W) 

Fig  13.  Half  -  tone  display  of  the  (a)  retrieved  cloud  temperature,  (b)  the  retrieved  ean  effective  size,  (c)  the  retrieved  optical  depth, 
and  (d)  the  retrieved  cloud  height,  over  the  5°  x  5°  area  over  Wisconsin  at  2100  UTC,  28  October  1986. 


The  computed  mean  effective  sizes,  however,  have  not  been  validated  against 
aircraft  observations,  because  the  aircraft  data  that  are  properly  collocated  and 
coincident  with  the  satellite  overpass  were  not  available.  In  Figure  13(c),  the 
retrieved  optical  depth  ranges  between  0  and  8 .  The  map  of  the  optical  depth 
matches  that  of  the  Ch.  4  brightness  temperature  shown  in  Figure  12.  In  Figure 
13(d),  the  distribution  of  cloud  heights  is  determined  from  the  distribution  of 
cloud  temperature  using  the  averages  of  the  2100  UTC  soundings  data  at  Plattville 
and  Green  Bay,  Wisconsin  (Starr  and  Wylie,  1990).  The  cloud  heights  range 
between  6  and  11  km.  The  distribution  of  cloud  height  agrees  with  that  of  cloud 
temperature  [Figure  13(a)]  and  mean  effective  ice  crystal  size  [Figure  13(b)]. 

The  preceding  results  are  compared  with  results  presented  by  Minnis  et  al . 
(1990b),  who  also  derived  the  distributions  of  cloud  height  and  cloud  optical 
depth,  using  both  the  AVHRR  visible  and  IR  channel  data  over  the  same  area  and 
for  the  same  time  period.  Minnis  et  al.  (1990b)  derived  the  optical  depth 
directly  from  the  Ch.  1  bidirectional  reflectances  based  on  radiative  transfer 
calculations  (Takano  and  Liou,  1989b).  The  present  removal -retrieval  scheme, 
however,  derives  the  optical  depth  based  on  the  parameterization  of  the  Ch.  4 
emissivity.  Minnis  et  al .  (1990b)  showed  that  the  optical  depth  ranges  between 
0  and  5  for  most  of  the  area.  Over  the  major  cloudy  areas  near  the  northwest 
corner  and  at  about  45®N  and  89®W,  the  optical  depth  is  shown  to  be  greater  than 
5.  The  distribution  of  the  optical  depth  from  the  present  removal -retrieval 
scheme  agrees  reasonably  well  with  that  presented  by  Minnis  et  al .  (1990b).  The 
present  results  show  more  detailed  features,  however.  The  distribution  of  cloud 
height  derived  by  Minnis  et  al.  (1990b)  is  between  7  and  11  km  for  most  areas. 
It  should  be  noted  that  the  cloud  heights  presented  by  Minnis  et  al .  (1990b)  are 
the  cloud-center  heights,  which  were  derived  through  empirical  formulas  that 
relate  cloud  temperature  and  cloud- top  temperature.  In  the  present  scheme, 
however,  the  retrieved  cloud  height  is  close  to  the  cloud- top  height  for  thicker 
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clouds.  Therefore,  over  most  of  the  cloudy  areas,  the  cloud  heights  determined 
from  the  present  scheme  are  slightly  higher  than  those  presented  by  Minnis  et  al. 
(1990b) . 

Finally,  for  verifications  and  comparisons,  Table  4  lists  the  retrieved 
values  of  cloud  temperature,  mean  effective  ice  crystal  size,  cloud  height,  and 
optical  depth  from  the  present  scheme.  Also  listed  are  the  cloud  temperature  and 
midcloud  height  derived  from  lidar  and  sounding  data,  and  the  optical  depth 
derived  from  the  GOES  data  (Minnis  et  al.  1990a).  Comparisons  are  made  near  the 
three  FIRE-I-IFO  stations  whose  locations  are  identified  in  Figure  12:  Madison 
(MAD,  43.1^N,  89.4®W),  Fort  McCoy  (FMC,  43.9®N,  90.8®W),  and  Wausau  (WAU,  45.0®N, 
89.7®W).  The  retrieved  values  from  the  present  scheme  are  based  on  the  average 
of  all  pixels  within  a  1°  x  1°  area  around  each  station.  Cloud  temperatures  and 
cloud  heights  at  all  three  stations  determined  from  the  present  scheme  agree 
reasonably  well  with  those  from  lidar  data  (Sassen  et  al .  ,  1990).  The 
differences  in  cloud  temperature  and  cloud  height  are  less  than  4.2  K  and  0.5  km, 
respectively.  These  differences  may  be  caused  by  the  effect  of  the  vertical 
inhomogeneity  within  cirrus  clouds.  The  optical  depths  inferred  from  the  present 
scheme  are  slightly  larger  than  those  derived  from  the  GOES  data  by  about  0.1. 
It  is  noted  that  Minnis  et  al .  (1990a)  derived  their  optical  depths  based  on  a 
cirrostratus  microphysical  model  with  a  mean  effective  size  of  about  41.5  ^ra. 
The  mean  effective  size  retrieved  from  the  present  study  is  between  50  and  60  jura. 
In  summary,  the  retrieved  cloud  parameters  agree  reasonably  well  with  the 
observations  and  with  results  computed  from  different  techniques. 

4.3  Retrievals  over  a  20°  x  30°  Area 

We  have  performed  the  daytime  retrieval  of  cirrus  cloud  parameters  covering 
the  midwestern  United  States  for  the  time  period  2100  UTC ,  28  October  1986. 
Figure  14  shows  the  display  of  Ch.  4  brightness  temperatures.  The  area  covered 
is  from  30°  to  50°N  and  from  80°  to  110°W,  It  is  evident  that  clouds  are  present 
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over  the  Great  Lakes  region  (outlined  by  dots  over  the  northeast  corner)  and  the 
southeast  United  States.  The  rest  of  the  area  is  relatively  clear.  For 

comparison,  Figure  15  displays  the  AVHRR  Ch.  1  satellite  imagery  based  on  HRPT 
(High-Resolution  Pixel  Transmission)  data  over  the  same  area  and  time  period. 

The  data  domain  is  divided  into  70  3®  x  3®  boxes.  For  each  box,  we  compute 
a  single  set  of  average  sun- satellite  geometrical  parameters  for  the  purpose  of 
correlating  Chs.  1  and  3  reflectances  based  on  radiative  transfer  calculations. 
The  variations  in  the  solar  zenith  angle,  satellite  zenith  angle,  and  relative 
azimuthal  angle  for  pixels  within  each  box  are  less  than  3®,  8®,  and  2®, 
respectively.  Radiative  transfer  calculations  show  that  these  angular  variations 
produce  negligible  variations  in  the  Chs.  1  and  3  reflectances.  Each  box  is 
further  divided  into  nine  1®  x  1®  mini -boxes.  The  clear/cloud  detection  program 
is  applied  to  all  pixels  in  each  mini -box  to  obtain  the  Ch.  1  effective  surface 
albedo  from  the  mode  of  the  histogram  of  all  Ch.  1  clear  reflectances  in  the  box. 
We  then  apply  the  removal -retrieval  scheme  to  cloudy  pixels  in  each  mini -box, 
utilizing  a  pair  of  Chs.  3  and  A  thermal  infrared  clear  radiances  for  each  mini¬ 
box,  If  no  clear  pixels  are  detected  within  the  mini-box,  the  clear  radiances 
of  a  neighboring  box  are  used.  The  averages  of  Green  Bay  and  Plattville 
soundings  at  2100  UTC  are  employed  to  compute  the  transmittances  for  clear 
atmospheres , 

Figures  16(a) -(c)  display  the  retrieved  cirrus  cloud  temperatures,  optical 
depths,  and  mean  effective  sizes,  respectively.  Over  the  northern  Great  Lakes 
region,  cloud  temperatures  are  lower  (-230  K)  ,  with  larger  pockets  of  cold  cirrus 
(<230  K)  than  over  the  southeastern  corner.  Moreover,  the  corresponding  cloud 
optical  depths  are  larger  than  8  over  the  Great  Lakes  region.  Small  mean 
effective  sizes  (<50  /urn)  are  found  near  the  northern  edge  of  Lake  Superior. 
These  retrieval  results  are  consistent  with  the  high-resolution  satellite  cloud 
picture  shown  in  Figure  15.  To  the  northeast  and  east  of  the  Great  Lakes,  the 
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Fig  1 4.  Display  of  the  Ch.  4  brightness  temperatures  over  the  midwestern  United  States  at  21 00  UTC,  28  October,  1 986. 


Fig  15.  The  AVHRR  Ch.  1  satellite  imagery  over  the  midwestern  United  States  at  2100 
UTC,  28  October  1986. 
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Fig  16.  Display  of  (a)  retrieved  cloud  temperature,  (b)  the  retrieved  optical  depth,  and  (c)  the  retrieved 
mean  effective  sizes,  over  the  area  shown  in  Fig.  8. 
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Table  4.  Cirrus  cloud  temperature,  mean  effective  ice  crystal  size,  cloud 
height,  and  optical  depth  determined  from  the  present  retrieval 
program. 


Tc(K) 

present  lidar* 

De(/im) 

present 

z(km) 

present  lidar* 

r 

present 

GOES* 

WAU 

226.5 

230.7 

57.1 

9.5 

9.0 

1.54 

1.43 

FMC 

229.7 

226.5 

61.3 

9.1 

9.5 

1.41 

MAD 

225.6 

228.1 

53.9 

9.6 

9.3 

0.6 

0.56 

*after  Minnis  et  al.  (1990a). 

area  is  mostly  covered  with  clouds.  The  visible  satellite  cloud  picture 
indicates  that  the  same  area  is  relatively  dark,  due  probably  to  the  large  solar 
zenith  angle  over  that  area  (-  80°).  Moreover,  the  Ch.  4  brightness  temperatures 
are  less  than  245  K  and  the  brightness  temperature  difference  between  Chs .  4  and 
5  are  between  2  and  3  K,  indicating  the  presence  of  thick  cirrus  over  that 
region.  Although  it  is  possible  that  low-level  clouds  could  have  been  present, 
reliable  schemes  for  detecting  multi-layer  cloud  systems  are  yet  to  be  developed 
to  identify  this  situation. 
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5. 


VERIFICATION  USING  FIRE-II-IFO  MEASUREMENTS 


5.1  Data  Sources 

The  FIRE-II-IFO  was  carried  out  near  Coffeyville,  Kansas,  during  November 
and  December,  1991.  One  of  the  objectives  of  this  field  experiment  was  to 
quantify  the  capabilities  and  limitations  of  methods  for  retrieving  physical  and 
optical  cirrus  cloud  properties  from  satellite  observations.  During  the  latter 
part  of  November  and  the  first  part  of  December,  1991,  a  total  of  11  cirrus  days 
was  observed.  For  all  these  cases,  we  have  reviewed  the  available  satellite  data 
and  lidar  measurements  as  well  as  examined  the  temporal  differences  between  the 
satellite  overpass  time  and  the  sampling  time  of  the  sounding  system.  We  find 
that  around  2100  UTC  on  both  26  November  and  5  December,  there  were  well  defined 
cirrus  cloud  layers  over  the  Coffeyville  area  and  that  the  temporal  differences 
between  the  satellite  overpasses  and  replicator  measurements  are  less  than  30 
minutes.  Thus,  we  have  focused  on  these  two  dates  for  the  retrieval  analyses  and 
verifications. 

For  the  retrieval  of  cirrus  cloud  properties,  high-resolution  (HRPT)  AVHRR 
data  from  NOAA-11  satellite  are  used.  This  satellite  has  nominal  equator 
crossing  times  of  0140  and  1340  LST.  The  AVHRR  is  comprised  of  five  channels: 
0.63,  0.8,  3.7  (NIR)  ,  10.9  (IR)  ,  and  12  /ira.  The  NIR  and  IR  radiances  were 
calculated  from  the  raw  counts  provided  in  the  NOAA-IB  data  stream  using  the 
nominal  calibration  (Kidwell  1991)  and  from  the  nonlinearity  corrections 
described  by  Weinreb  et  al .  (1990).  The  visible  channel  radiances  were 
calculated  using  the  instrument  degradation  correction  function  provided  by 
Whitlock  et  al .  (1990). 

Balloon-borne  formvar  (polyvinyl  formal)  ice  crystal  replicators  were 
launched  during  FIRE-II-IFO  to  measure  the  "vertical  profiles"  (in  a  Lagrangian 
sense)  of  cirrus  microphysical  properties,  with  emphasis  on  the  detection  of 
small  ice  particles.  Schaefer  (1941)  first  developed  the  formvar  replication 
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technique  to  capture  and  preserve  detailed  features  of  snowflakes,  using  a 
solution  of  formvar  dissolved  in  ethylene  dichloride.  Particles  falling  onto  a 
solution-coated  slide  or  filmstrip  are  enveloped  by  fluid,  forming  a  detailed 
replica  of  the  particle  when  the  solvent  evaporates  and  the  plastic  hardens.  The 
use  of  formvar  replication  on  a  balloon-borne  instrument,  i.e.  a  "replicator 
sonde”,  was  developed  by  Magono  and  Tazawa  (1966)  to  study  the  vertical 
distribution  of  snow  crystals  up  to  500  mb  level  in  a  snow  cloud.  The 
fundamental  design  of  the  ice  crystal  replicator  employed  in  FIRE-II-IFO  was 
based  on  the  prototype  developed  by  Magono  and  Tazawa,  with  modifications  to 
improve  the  data  quality  and  to  allow  the  instrument  to  collect  data  in  high  and 
cold  cirrus  clouds. 

Replicator  launches  were  timed  to  coincide  with  satellite  overpasses  and 
during  the  time  periods  when  ground-based  remote  sensors  and  aircraft  were 
operating.  The  replicators  were  launched  together  with  the  NCAR- CLASS  and  then 
tracked  and  retrieved  by  an  automobile  using  the  radio  direction- finding  (RDF) 
equipment.  The  balloon  size  and  inflation  were  chosen  to  deliver  an  ascent  rate 
of  about  5  m/sec  and  to  burst  at  11-13  km  altitude,  after  which  the  train 
descended  on  a  parachute.  More  comprehensive  discussions  on  the  replicator  data 
have  been  presented  by  Heymsfield  and  Miloshevich  (1993).  High  resolution 
radiosonde  data  and  position  coordinates  were  obtained  from  the  NCAR-CLASS. 

5.2  Results  and  Discussions 

As  described  previously,  the  Fire-II-IFO  over  Coffeyville  was  carried  out 
during  November -December  1991  with  a  total  number  of  11  observed  cirrus  days 
(November  13,  17,  22,  25,  26,  27,  28,  29,  30,  and  December  5  and  6).  Based  on 
inspection  of  the  available  AVHRR  data,  replicator  measurements,  NCAR-CLASS 
sounding  data,  and  surface  broadband  radiative  flux  measurements,  we  find  that 
the  data  collected  during  5  December  are  most  appropriate  for  the  present 
retrieval  analyses  and  verifications  for  the  following  reasons.  First,  on  that 
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date,  a  layer  of  well  defined  cirrus  cloud  was  present  over  Coffeyville,  so  that 
the  retrieval  scheme,  which  was  principally  designed  to  retrieve  the  single- layer 
cloud  properties,  can  be  directly  applied  to  the  satellite  data.  Second,  the 
temporal  differences  between  the  satellite  overpass  time  and  the  replicator 
sampling  time  are  small.  The  two  data  sets  are  therefore  closely  coincident  and 
can  be  collected  in  a  Lagrangian  sense.  This  part  will  be  discussed  further  for 
each  case.  Finally,  the  satellite  viewing  angles  are  also  small  (<  30°).  This 
condition  guarantees  that  the  size  for  each  pixel  is  small,  so  that  the  effects 
of  cloud  horizontal  inhomogeneity  on  the  analysis  can  be  minimized. 

On  5  December,  the  sky  was  clear  at  sunrise  and  not  very  cold.  A  rapid 
influx  of  light  cirrus  (Spissatus)  was  observed  from  the  west  during  the  morning. 
Coverage  and  density  of  cirrus  increased  through  the  day.  Aircraft  and  ground- 
based  sensors  estimated  cloud  bases  in  the  10  km  range,  with  tops  extending  to 
around  13  km.  During  the  late  afternoon,  multi-layered  features  covered  the  sky, 
with  all  stations  in  Kansas  reporting  high  level  clouds.  No  low  or  middle  level 
clouds  passed  over  Coffeyville,  but  some  low  clouds  were  visible  in  the  distant 
southern  sky  during  the  afternoon. 

Figures  17 (a) -(c)  show  the  distribution  of  Ch.  4  brightness  temperatures, 
and  the  retrieved  optical  depths  and  mean  effective  ice  crystal  sizes  over  the 
0.5°  X  1.0°  domain  (36.75°  -  37.25°  N;  94.8°  -  95.8°  W)  at  2100  UTC,  5  December 
1991,  respectively.  The  symbol  denotes  the  location  of  Coffeyville.  It  is 
noted  that  almost  the  whole  domain  was  covered  with  cirrus  clouds.  The 
distribution  of  the  retrieved  optical  depths  closely  match  that  of  the  Ch.  4 
brightness  temperature.  The  cloud  optical  depth  northeast  of  Coffeyville  and 
near  the  southern  boundary  of  the  domain  is  larger  than  3.  However,  the 
distribution  of  mean  effective  ice  crystal  sizes  show  no  correlation  with  that 
of  the  Ch.  4  brightness  temperature.  Generally,  the  mean  effective  size 
increases  from  the  south  toward  the  north.  The  reliability  of  these  retrieval 
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Fig  17.  Display  of  (a)  AVHRR  Ch.  4  brighness  temperature,  (b)  the  retrieved  optical  depth,  and  (c)  the 
retrieved  mean  effective  ice  crystal  size  over  a  0.5®  x  1 .0®  area  around  Coffeyville,  Kansas  (denoted 
by  the  symbol  V),  at  2108  UTC,  5  December  1991 . 


53 


results  can  only  be  confirmed  by  comparing  with  independently  derived 
observational  data . 

Figures  18 (a) -(e)  show  the  ice  crystal  size  distributions  derived  from  the 
replicator  data  at  five  selected  height  levels  between  the  cloud  base  and  cloud 
top.  Near  the  cloud  top  (~  13  km),  the  cirrus  cloud  is  composed  of  small  quasi- 
spherical  particles  (L  <  100  /im) .  The  smallest  ice  crystals  that  the  replicator 
can  measure  are  on  the  order  of  5  /zm.  In  the  middle  part  of  the  cloud  (9.8  - 
12.0  km),  there  are  columns,  bullets,  plates,  and  irregular  particles.  The  size 
distributions  become  widen  and  more  large  particles  appear  in  the  lower  part  of 
the  cloud.  Near  the  9. 77 -km  level,  the  size  distribution  is  bimodal,  where  the 
larger  mode  is  presumably  due  to  the  presence  of  large  particles  falling  out  from 
higher  levels.  Finally,  near  the  cloud  base  (-9.6  km),  there  are  only  partially 
sublimated,  irregularly  shaped  ice  crystals. 

Because  clouds  are  horizontally  inhomogeneous,  collocation  and  coincidence 
of  the  satellite  data  and  replicator  measurements  in  a  Lagrangian  sense  are 
required  for  validation  of  the  retrieval  results.  Figure  19  shows  the 
geographical  coordinates  of  the  replicator  at  various  times  after  launching.  The 
replicator  was  launched  at  2045  UTC  and  reached  the  cloud  base  and  top  at  2113 
and  2125  UTC,  or  5  and  17  minutes  after  the  satellite  overpass,  respectively. 
These  temporal  differences  are  relatively  small  compared  to  the  life  time  of 
cirrus  clouds,  which  is  on  the  order  of  several  hours.  To  collocate  the  two  data 
sets,  we  must  determine  the  precise  position  of  the  cloudy  region  under  the 
satellite  overpass  that  was  sampled  by  the  replicator.  In  Figure  19,  the  open 
and  closed  circles  denote  the  coordinates  of  the  sounding  system  when  it  reached 
the  cloud  base  and  top,  respectively.  Using  the  cloud- level  windspeed  and  wind 
direction  data  recorded  by  the  NCAR- CLASS  sounding  system,  we  determine  that  at 
the  time  of  satellite  overpass  the  cloud  associated  with  the  replicator 
measurements  was  located  at  a  position  denoted  by  "x" .  This  point 
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Fig  18.  Ice  crystal  size  distributions  derived  from  the  replicator  data  at  five  selected  height 
levels  within  the  cloud,  which  was  sampled  between  2113  and  2125  UTC,  5  December 
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Fig  1 9.  Geographical  coordinates  of  the  replicator  at  various  times  after  launching  on  5  December  1 991 . 


near  Coffeyville  was  upwind  of  the  sounder  positions  at  the  cloud  base  and  top, 
because  the  satellite  overpass  was  before  the  replicator  sampling.  Since  the 
coordinates  of  this  point  are  subject  to  uncertainties  in  the  wind  data,  an  area 
of  0.05*^  X  0.2°  around  this  point  is  subsequently  selected  for  the  verification 
s  tudy . 

Figures  20(a) -(b)  show  the  averaged  retrieval  results  within  the  0.05°  x 
0.2°  domain  along  with  sounding  measurements.  Figure  20(a)  displays  the 
temperature  and  relative  humidity  profiles  obtained  from  the  NCAR-CIASS  sounding 
system  between  2045  and  2129  UTC.  The  cloud  base  and  top  heights  derived  from 
the  replicator  sounding  are  9.5  and  13  km,  respectively.  The  relative  humidity 
sharply  increases  with  height  near  the  cloud  base  and  decreases  upward  above  10 
km.  The  mean  retrieved  cloud  temperature  is  212  K,  which  is  the  average  of  81 
pixels  within  the  0.05°  x  0.2°  area.  The  standard  deviation  is  3  K,  revealing 
that  cloud  temperatures  are  relatively  uniform  within  the  retrieval  domain. 
Moreover,  the  mean  cloud  height  determined  from  the  temperature  sounding  is  11.5 
km,  which  is  within  the  observed  cloud  boundaries.  It  should  be  noted  that  the 
mean  Ch.  4  brightness  temperature  is  255  K,  much  higher  than  the  mean  retrieved 
cloud  temperature.  If  this  temperature  were  used  to  retrieve  the  cloud  height, 
the  error  will  be  on  the  order  of  6  km. 

Figure  20(b)  shows  the  vertical  distribution  of  mean  effective  ice  crystal 
sizes  computed  from  the  replicator  data.  The  mean  effective  size  is  15  /im  near 
the  cloud  top  and  increases  to  135  near  the  cloud  base.  From  these  values, 
we  determine  the  vertically  averaged  mean  effective  ice  crystal  size  to  be  97.3 
/zm,  as  shown  in  Figure  20(b)  by  the  solid  vertical  bar.  This  size  is  the  sum  of 
the  mean  effective  ice  crystal  size  at  a  given  level  weighed  by  the  number 
density  at  that  level  normalized  by  the  total  number  density.  The  retrieved  mean 
effective  size  is  88.7  /zm,  which  differs  from  the  replicator  value  by  less  than 
10  /im.  On  the  bottom  scale  of  Figure  20(b)  are  shown  the  replicator-derived  and 
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the  retrieved  optical  depths,  which  are  1.38  and  1.49,  respectively.  The 
differences  in  the  retrieved  and  replicator-derived  mean  effective  size  and 
optical  depth  could  well  be  due  to  the  small  time  difference  between  satellite 
pass  and  replicator  sounding. 
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6. 


CONCLUSIONS 


A  removal -retrieval  technique  based  on  radiative  transfer  parameterizations 
has  been  developed  for  inferring  cirrus  cloud  parameters  during  local  daytime. 
This  scheme  uses  radiance  data  of  AVHRR  3.7  and  10.9  urn  channels  to 
simultaneously  determine  cirrus  temperature,  mean  effective  ice  crystal  size,  and 
optical  depth.  During  the  local  daytime,  the  3.7  /xra  (Ch.  3)  solar  component  is 
removed  by  using  the  correlation  between  Ch.  1  (0.63  ^m)  and  Ch.  3  solar 
reflectances.  This  correlation  is  derived  from  radiative  transfer  calculations. 

We  have  carried  out  sensitivity  studies  on  the  accuracy  of  the  removal 
scheme.  Overall,  for  optical  depths  larger  than  1,  the  maximum  percentage  error 
in  the  3.7  pm  solar  component  is  less  than  10%.  Moreover,  we  have  performed 
retrievals  over  a  1°  x  1°  area  west  of  Fort  McCoy,  Wisconsin,  using  the  AVHRR 
data  collected  at  2100  UTC,  28  October  1986.  The  distribution  of  Ch.  4 
brightness  temperature  over  this  area  shows  distinct  clear  regions.  Based  on  the 
retrieved  cloud  temperature,  several  cold  cirrus-generating  regions  are  shown 
along  with  large  areas  of  warm,  stable,  and  dissipating  cirrus  clouds  that  are 
present  to  the  north  and  east  of  the  cirrus  -  forming  regions.  Significant 
differences  are  evident  for  the  patterns  of  the  retrieved  optical  depth  and  cloud 
temperature.  The  retrieved  mean  effective  ice  crystal  sizes  range  from  30  to  90 
pm,  which  appear  to  be  consistent  with  the  observed  values. 

We  have  applied  the  removal-retrieval  scheme  to  retrieve  cirrus  parameters 
over  a  5®  X  5®  area,  using  the  same  daytime  data  set  collected  on  28  October 
1986.  The  retrieved  cloud  temperatures  and  cloud  heights  compare  reasonably  well 
with  the  results  presented  by  other  researchers  using  different  methods  and  with 
available  lidar  and  aircraft  measurements.  We  have  also  carried  out  the  mapping 
of  cirrus  cloud  parameters  over  the  midwestern  United  States .  The  distributions 
of  the  retrieved  cirrus  cloud  parameters  are  consistent  with  the  satellite  cloud 
picture  based  on  the  high-resolution  AVHRR  Ch.  1  reflectances. 
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Moreover,  verifications  of  the  retrieved  cirrus  optical  depth  and  ice 
crystal  sizes  have  been  performed  using  balloon-borne  replicator  and  sounding 
data  obtained  during  FIRE-II-IFO.  We  have  taken  the  advantage  that  the  balloon- 
borne  replicator  data  can  provide  a  nearly  continuous  vertical  record  of  ice 
crystal  size  distributions  in  a  Lagrangian  sense.  An  analytical  method  has  been 
developed  to  derive  the  optical  depth  and  the  mean  effective  ice  crystal  size 
from  the  replicator  data,  taking  into  account  the  effects  of  different  shapes, 
including  quasi-spherical ,  irregular,  and  hexagonal  ice  crystals.  Moreover,  we 
apply  the  retrieved  cloud  properties  to  the  computation  of  surface  radiative 
fluxes  using  a  radiative  transfer  program  that  involves  a  consistent 
representation  of  cloud  fields  and  compare  the  computed  values  with  data  measured 
from  ground-based  radiometers . 

We  have  focused  on  the  case  of  5  December.  On  that  date,  a  well  defined 
cirrus  layer  was  present  over  the  Coffeyville  area  and  the  temporal  differences 
between  the  satellite  overpass  and  the  replicator  sampling  are  small.  The  two 
data  sets  are  closely  coincident  and  can  be  collocated  in  a  Lagrangian  sense. 
In  each  case,  we  first  perform  retrievals  over  a  0.5°  x  1.0°  region  around  the 
Coffe3rville  area  and  obtain  the  distribution  of  cloud  optical  depth  and  mean 
effective  size.  We  then  determine  the  precise  position  of  the  cloudy  region  that 
was  sampled  by  the  replicator  under  the  satellite  overpass.  A  small  area  around 
this  position  is  subsequently  chosen  and  the  mean  retrieved  cloud  parameters  of 
this  domain  are  compared  with  those  derived  from  the  replicator  and  sounding 
data.  We  show  that  the  retrieved  cirrus  cloud  temperature,  mean  effective  size, 
and  optical  depth  closely  match  with  the  values  determined  from  sounding  and 
replicator  data.  Differences  in  the  retrieved  and  replicator- derived  values  of 
the  cloud  parameters  could  be  a  result  of  the  small  temporal -differences  between 
the  satellite  pass  and  replicator  sounding. 
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